Optimal broadening of finite energy spectra in the numerical renormalization group: 
application to dissipative dynamics in two-level systems 
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Numerical renormalization group (NRG) calculations of quantum impurity models, based on a 
logarithmic discretization in energy of electronic or bosonic Hamiltonians, provide a powerful tool to 
describe physics involving widely separated energy scales, as typically encountered in nanostructures 
and strongly correlated materials. This main advantage of the NRG was however considered a 
drawback for resolving sharp spectral features at finite energy, such as dissipative atomic peaks. 
Surprisingly, we find a bunching of many-body levels in NRG spectra near dissipative resonances, 
and exploit this by combining the widely-used Oliveira's z-trick, using an averaging over few discrete 
NRG spectra, with an optimized frequency-dependent broadening parameter b{uj). This strategy 
offers a tremendous gain in computational power and extracts all the needed information from the 
raw NRG data without a priori knowledge of the various energy scales at play. As an application 
we investigate with high precision the crossover from coherent to incoherent dynamics in the spin 
boson model. 
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A general hallmark of many-particle interaction, as 
found in a variety of condensed matter systems such 
as nanostructures and strongly correlated materials, lies 
in the presence of several energy scales, possibly widely 
separated from each other due to renormalization and 
dynamical effects. Two well-studied examples found in 
electronic systems concern the Kondo effect for magnetic 
impurities in metals and Fermi liquids in proximity to a 
Mott insulating phase, two instances where low-energy 
quasiparticles emerge below a typical temperature which 
is quite reduced from the bare Fermi energy^. These low- 
lying excitations do however coexist with higher energy 
atomic levels, also broadened and displaced in a strong 
manner from their bare atomistic values due to the dis- 
sipation brought by the electronic environment. Such 
complex physical effects, taking place on a broad range of 
energies, entail great practical difficulties for most direct 
numerical approaches. These are partially lifted using 
Wilson's original idea of the logarithmic discretization3i2,, 
as implemented in numerical renormalization group cal- 
culations (see Rcf. Q for a recent review) . This technique 
has been improved in the last twenty years to calculate 
static and dynamic quantities both for fermionio^ and 
bosonic quantum impurity models^. Important practical 
applications until now involve the calculation of trans- 
port in the Kondo regime for Kondo alloys and artifi- 
cial quantum dots^, as well as the accurate description 
of the zero-temperature Mott transition by combining 
NRG^ with Dynamical Mean Field Theory (DMFT)S. 
More generally, exponentially small energy scales are also 
found in the vicinity of quantum critical points^^, so that 
impurity models provide a simplified testbed for the the- 
ory of quantum critical phcnomenaii. Again the NRG 
is the ideal technique for studying such impurity quan- 
tum phase transitions^, with potential implications for 
artificial nanostructuresi^ii^. 

Despite these successes, the foundation of NRG on a 



logarithmic discretization mesh implies a loss in accuracy 
for the high energy spectral features, which has plagued 
most calculations so far. Not only are atomic-like excita- 
tions physically observable, using photoemission or tun- 
neling spectroscopies, but they may also be intimately 
tied to the low-energy excitations. Such interesting be- 
havior occurs in the vicinity of the Mott metal-insulator 
transition?., where the self-consistcntly determined elec- 
tronic environment in the DMFT picture shows electronic 
states violently redistributed over all energy scales. The 
numerical cost of converging the DMFT equations cer- 
tainly requires efficient and accurate NRG calculations 
for the spectral functions, without a priori knowledge of 
the excitations involved. For this reason, the idea of av- 
eraging over Nz realization of the Wilson chain, using 
the so-called z-trick^^ to fill in the missing spectral infor- 
mation (to be discussed later on), seems prohibitive for 
most practical calculations, and offers a limited gain for 
very narrow spectral structures (see Ref. [l^ for a detailed 
study). 

In this Rapid Communication, we show that the broad- 
ening procedure used to smoothen the discrete NRG data 
is one critical step for obtaining an optimal resolution at 
all energies. We henceforth develop a simple adaptive 
procedure where the standard broadening parameter h 
is taken to be frequency dependent. This choice is dic- 
tated by our surprising observation that the density of 
many-body NRG levels increases sharply as soon as nar- 
row atomic resonances are encountered. Together with 
the usual z-trick, combining several but reasonably few 
NRG calculations, this extra &(cj)-trick allows to com- 
pute finite frequency spectra using as little as = 10 
NRG calculations, in situations where large scale NRG 
z-averaging would be prohibitive. An improved broad- 
ening procedure leading to errors in spectral functions 
limited to few percents could constitute a further step in 
applying the NRG to a wider class of problems, such as 
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DMFT+NRG calculations. 

In order to explicitly demonstrate these ideas, we focus 
on the simplest quantum impurity model, namely the 
spin boson Hamiltonian: 



(1) 



that involves a two-level system, described by a quantum 
spin 1/2, and a bosonic bath a\ with continuous spectrum 
LOi of energies. A is a transverse magnetic field, while the 
coupling constant A controls the strength of longitudinal 
dissipation. The bosonic spectrum is assumed to be sub- 
ohmic with bath exponent < s < 1 (this includes the 
well-studied ohmic case s = 1): 

J{uj) = ^Ti\^5{uo-uj,) = 2Tiaujl~''uo'e{uo)e{uOc-i^) (2) 



where lOc is a high-energy cutoff. For small values of the 
dissipation a, this model is known to exhibit coherent 
precession of the spin around the x axis at zero tempera- 
ture. By increasing the coupling to the bath for values of 
s < 1, Rabi oscillations arc progressively damped^^, be- 
fore the two-level system localizes in one potential min- 
ima via a quantum phase transition^. When s ^ how- 
ever, the phenomena of localization and decoherence oc- 
cur in reverse order—. Despite its simplicity, this model 
embodies all the effects typical of strong correlations: 
low-energy scales are indeed dynamically generated near 
the localization transition, while discrete spin excitations 
due to the transverse field A are deeply affected by the 
bosonic environment, that leads to broadening and fre- 
quency shift in the magnetic response 
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(3) 



with C{t) = \{[(T^{t),a^]) the spin autocorrelation func- 
tion. 

The implementation of the NRG procedure follows the 
standard routed as discussed in Ref. @ for the exten- 
sion to bosonic Hamiltonians. The bosonic band ^ 
is logarithmically discretized with the Wilson parameter 
A > 1, first on the highest energy interval [K~^ijJc,iOc\, 
and then iteratively on successive decreasing energy win- 



dows [A 



c] for n strictly positive in- 



teger. This also introduces the crucial z parameter 
(0 < z < 1) that is used to average over Wilson 
chainsi^, allowing to obtain better resolution on the fi- 
nite energy states. The rest of the NRG follows Rcf. iQ, 
coupling iteratively the kept energy levels up to iteration 
n to the states living in the shell n + \, and truncating 
the successive Hamiltonians to keep up with a manage- 
able number of eigenstatcs. All subsequent calculations 
were performed with A = 2, A^j, = 8 kept bosonic states 
on the added bosonic "site", and A^ = 160 kept NRG 
states, ensuring good convergence. The resulting discrete 
spectra at successive NRG iterations are combined using 



the interpolation scheme proposed in Ref. [T^ (see how- 
ever Ref. hj for a more rigorous implementation) , leading 
to a set of z-dependent many-body energy levels ea,z la- 
belled by quantum number a. The spin-spin correlation 
function is thus readily obtained at zero temperature as: 
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where eo,2 is the ground state energy, and these raw NRG 
data arc displayed on Fig. [U to be discussed below. 
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FIG. 1: (Color online) Spin susceptibility C(lj) of the sub- 
ohmic spin boson model at s = 0.1, A = Q.lujc and oi = 
0.000125. Raw data \{Q,z\a^\a,z)\^ /[2N^{ea.z - eo,z)] are 
given as circles for Nz = 20 combined NRG calculations, solid 
line is the perturbative result ([6]), and the three dashed lines 
are the various NRG broadenings discussed in the text. 

Let us now discuss the general structure of the NRG 
spectra. For a single NRG calculation with a given value 
of z, the a priori energy resolution at the scale l^jn = 
A-^-^ is 6ujN = A-^-^ - A-^-^-i = (1 - A^ijwAT, 
so that the resolution degrades at increasing energy. For 
obtaining smooth NRG spectra, the delta-peaks in Q 
are usually broadened^ at energy ujpf on the same scale: 



S (\uj\ - un) 



-[log(|t^|/t^jv)V&] 



(5) 



with 6 w 0.7 typically. Combining NRG runs us- 
ing the z-averaging in ([T]) allows in principle to improve 
by a factor the accuracy at high energy, since the 
broadening parameter may now be decreased down to 
b = 0.7 /Nz- This procedure actually faces two problems: 
i) for very sharp resonances (typically several orders of 
magnitude narrower that the natural high-energy cut- 
off), parallelizing ^ 10 NRG calculations becomes 
too prohibitive, especially with the aim of DMFT sim- 
ulations; ii) states far from the resonances can then be 
too much under broadened, so that oscillations of period 
A/Nz can become quite prominent. In this view, quanti- 
tative spectra can be accurately extracted from the NRG 
data only in the continuum limit, either with A ^ 1 or 
Nz +00, as assumed in the previous litcraturei^i^iSi. 
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FIG. 2: (Color online) Comparison of various broadening pro- 
cedures near the resonance for the parameters of Fig. [T] stan- 
dard interleaved averaging with = 20 and b = 0.7 /N^ — 
0.035 (dash-dotted line), fe(a;)-trick with bo = 0.12 (dashed 
line), and analytical result (full line). The adaptive broaden- 
ing parameter b{ui) is also given (note the logarithmic scale 
on the right). 



We wish to show however that NRG provides much 
more information close to dissipative resonances than 
usually expected. For this purpose, we closely exam- 
ine the raw NRG spectra for resonance widths as small 
as 10^'^cjc, using reasonably few interleaved averaging 
Nz = 20, much smaller that the naively needed 10'' NRG 
runs. Fig. [T] shows that for frequencies far from the reso- 
nance located at A, the highest weight NRG states come 
in packets of Nz peaks for a given Wilson shell. Sur- 
prisingly, the density of NRG eigenstates demonstrates a 
huge increase precisely at the peak value, so that informa- 
tion about the resonance width and height seems really 
encoded in the discrete results! The standard broadening 
b = 0.1 /Nz = 0.035 is however way too large to benefit 
from this effect, and the corresponding smoothed spectra 
are indeed completely inaccurate. For definite compari- 
son, we have plotted the analytical result obtained from 
the expansion at weak dissipation^^.: 
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In order to judiciously exploit this unexpected finding, 
the straightforward idea we follow here is to adapt the 
broadening parameter 6 in ([5]) to the frequency depen- 
dence of the local density of NRG peaks. A very natural 
approach adapted to the NRG logarithmic discretization 
is to extract h{u)) from the logarithmic derivative of the 
integrated spectrum up to frequency lo: 



dlogw 



dlogT C 
dlogu; 



(7) 

where g ~ 1 is regularization parameter, whose precise 
value is not sensitive to the final results, and bo sets the 



typical broadening at very low and very high frequen- 
cies (far from the atomic resonances). Note indeed that 
we average in this expression two frequency sweeps from 
uj = and to = +oo respectively, in order to treat on 
an equal basis low and high frequency tails. Because the 
actual NRG data is fully discrete ([3]), we extract b{u!) re- 
cursively using equation ([7]) on the broadened NRG spec- 
tra. This procedure converges after few iterations to the 
results displayed on Figs. [1] and [2l Given the small nu- 
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FIG. 3: (Color online) Similar comparison as done on Fig. [2] 
with stronger dissipation a — 0.02. Optimized broadenings 
with 6(ii;)-trick were done with A^, — 2,20 and bo ~ 0.6,0.12 
respectively. 

merical effort devoted to generate these NRG data, the 
quality of final result is quite astonishing. The reason 
for this success is given on Fig. [21 where the b{uj) pa- 
rameter takes values as small as 10"'' at the resonance, 
while it increases drastically far away from it, naturally 
cancelling out a great part of the NRG oscillations due 
to the discretization of the Wilson chain. We emphasize 
that the only free parameter is the typical low-frequency 
broadening bo in which is easily adjusted: we al- 
ways find a range of bo values where the self-consistently 
broadened spectrum is relatively independent of the cho- 
sen bo. More generally, one can also avoid the use of 
Ansatz (O by determining b{uj) from estimates of the lo- 
cal density of highest weight states in the raw NRG data, 
giving similar results. Large scale NRG simulations with 
Nz ^ 10 and fixed b, as performed e.g. ini^, may thus 
benefit greatly from an adaptive broadening procedure. 
We also note that the high frequency tails are also faith- 
fully reproduced using the 6(u;)-trick. Finally, in situa- 
tions where relatively broader resonances are present, we 
show that quantitative results can be obtained at lower 
computational cost, by decreasing significantly the num- 
ber of z-averaging even as low as Nz < 5, see Fig. [31 

In order to check that our broadening procedure is 
completely robust for the whole range of parameters in 
the model (P), we investigate the effect of increasing dis- 
sipation for several bath exponent values s. The quanti- 
tative comparison of the Rabi resonance to the analytical 
formula ^ in its domain of validity, namely a <C 1 , gives 
us indeed confidence in the adaptative method. This is 
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proved on Fig. ffland [5l which consider the same param- 
eters as in Ref. llTl . with important quantitative improve- 
ment. For the smah value of s = 0.1 taken on Fig.lH the 



non-perturbative effects on the resonance are seen inbe- 
twecn the weak couphng regime and the quantum critical 
point, so that perturbation theory in a breaks down. 
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FIG. 4: (Color online) Comparison of optimized NRG spec- 
tra (solid lines) and perturbative calculations (dashed lines) 
with bath exponent s — 0.1, A = O.luJc and increasing values 
(arrow) of dissipation a = 0.002, 0.004, 0.007. 




FIG. 5: (Color online) Similar plot as on Fig. |4l with bath 
exponent s = 0.5, A = O.lui^ and a = 0.02, 0.06, 0.10. 



perfect matching of the NRG data on the resonance to the 
lowest order calculation in a shows that the dissipation 
mechanism does not care for the low-energy behavior of 
the spin dynamics, even at the quantum critical point^iii 
where the spin localizes for the values etc ~ 0.0071 and 
0.105 with s = 0.1 and 0.5 respectively {A/uc — 0.1 and 
A = 2 here) . Interestingly, the low-energy tails show sig- 
nificant deviations from the analytical result ^ even far 
from the quantum phase transition, so that lowest-order 
calculations^^ do not apply for the long-time dynamics. 
Improved resummation of the perturbation theory to all 
orders in a will be considered in a forthcoming workS^. 
For the intermediate value s = 0.5 shown on Fig.[5l clear 



To conclude, we have investigated the properties of dis- 
crete NRG spectra near atomic states in the spin boson 
model, and found that a bunching of the many-body lev- 
els occurs whenever sharp resonances are encountered. 
An adaptive broadening scheme was proposed, showing 
a drastic improvement in computation power for the cal- 
culation of accurate spectral functions over the whole en- 
ergy range. This procedure will certainly allow to take 
further advantage of the potentialities of the NRG in a 
wide range of physical situations, from ab-initio calcula- 
tions for magnetic impurities in metals^l to the difficult 
problem of simulating strongly correlated materials in the 
framework of the Dynamical Mean Field Theory^. 
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